function [out_struct] = combine_mat(information,syn_sim,kinematics)

% subject information
out_struct.subject = information.subject;
out_struct.eva = information.eva;
out_struct.task = information.task;

% kinematics
out_struct.pkvel = kinematics.hand_peak_velocity_mean;
%out_struct.handvel_fit_R2 = kinematics.handvel_gaussian_R2_mean ;
out_struct.bell_error = kinematics.bell_error_mean;
out_struct.hand_endpoint = kinematics.hand_endpoint_std;
out_struct.duration_time = kinematics.duration_time_mean ;
% out_struct.traj_sim_cos = kinematics.traj_sim_Cos_mean;
% out_struct.traj_sim_RMSe = kinematics.traj_sim_RMSe_mean;
out_struct.shoulder_flex = kinematics.shoulder_flextion_angle_mean;
out_struct.elbow_flex = kinematics.elbow_flextion_angle_mean;
% out_struct.acromine_excursion = kinematics.acromine_endpoint_mean;

% synergy sim
for i = 1:1:length(syn_sim)
    sim_muscle = syn_sim{i,1}.muscle_vector;
    sim_time = syn_sim{i,1}.time_profile;
    sim_combine = syn_sim{i,1}.syn_combine;
    switch i
        case 1
            out_struct.sim2_vector = sim_muscle;
            out_struct.sim2_timeprofile = sim_time;
            out_struct.sim2_combine = sim_combine;
        case 2
            out_struct.sim3_vector = sim_muscle;
            out_struct.sim3_timeprofile = sim_time;
            out_struct.sim3_combine = sim_combine;
         case 3
            out_struct.sim4_vector = sim_muscle;
            out_struct.sim4_timeprofile = sim_time;
            out_struct.sim4_combine = sim_combine;
         case 4
            out_struct.sim5_vector = sim_muscle;
            out_struct.sim5_timeprofile = sim_time;
            out_struct.sim5_combine = sim_combine;
    end
end

